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ABSTRACT 

Dynamical studies of superbubbles and Wolf-Rayet ring nebulae show discrepancies from the 
standard, adiabatic model for wind-blown bubbles. We therefore study the physical properties 
and kinematics of three candidate bubbles blown by single stars, to evaluate whether these 
discrepancies are also found in these simpler objects. Our sample candidates are N44F, N44J, 
and N44M, in the outskirts of the Hn complex N44 in the Large Magellanic Cloud. We have 
obtained ground-based and HST emission-line images and high dispersion echelle spectra for 
these objects. From the Ha luminosities and the [O in] /Ha ratios of these nebulae, we estimate 
the spectral types of the ionizing stars to be 07V, 09.5V and 09.5V for N44F, N44J, and 
N44M, respectively. We find that the observed expansion velocity of 12 km s _1 for N44F is 
consistent with the stellar wind luminosity expected from the central ionizing star, as predicted 
by the standard bubble model. The observed upper limits for the expansion velocities of N44 J 
and N44M are also compatible with the expected values, within the uncertainties. We also 
report the discovery in N44 F of strongly-defined dust columns, similar to those seen in the Eagle 
Nebula. The photoevaporation of these dense dust features may be kinematically important and 
may actually govern the evolution of the shell. The inclusion of photoevaporation processes may 
thus undermine the apparent agreement between the observed bubble dynamics and the simple 
adiabatic models. 
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Mechanical feedback from massive stars is a 
fundamental driver of galaxy evolution. The su- 
pernovae (SNe) and supersonic stellar winds from 
these stars generate bubbles and shells in the in- 
terstellar medium (ISM), which may dominate 
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the ISM structure formation in many star-forming 
galaxies. The standard, adiabatic model for these 
shells and superbubbles (e.g., Weaver et al. 1977) 
predicts the production of hot, 10 6 — 10 7 K, low- 
density gas within these shells, which bear the 
products of SN and massive star nucleosynthesis. 
Therefore, these bubble structures are thought to 
be the source of the diffuse, hot ionized compo- 
nent of the ISM. Vigorous star-forming regions 
may generate supergiant shells that blow out of 
galaxy disks, thereby dispersing metals and mass 
far from the parent star formation event. In ad- 
dition, this mechanical energy contributes to ISM 
kinematics and is thought to be a primary source 
of turbulence (e.g., Norman & Ferrara 1996; Gold- 
man 2000). 

The effects of mechanical feedback have been 
clearly demonstrated, for example, in starburst 
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galaxies (Johnson et al. 2000; Strickland et al. 
2000) and in superbubbles around OB associa- 
tions (e.g., Saken et al. 1992; Oey 1996). How- 
ever, several troubling discrepancies remain un- 
solved. For example, almost all clean examples 
of bubbles blown by Wolf-Rayet stars and of su- 
perbubbles present shells that are too small for 
their observed parent stellar population (e.g., Tr- 
effers & Chu 1982; Drissen et al. 1995). Many su- 
perbubbles also show anomalous kinematics and 
X-ray emission (e.g., Rosado et al. 1982; Chu & 
Mac Low 1990; Wang & Helfand 1991). To resolve 
these puzzles and more clearly understand the me- 
chanical feedback process, it is therefore necessary 
to examine simpler systems. 

In active star-forming regions, high concentra- 
tions of OB stars collectively produce superbub- 
bles. But in their peripheries, where stars are 
loosely distributed, discrete bubbles may be pro- 
duced by individual stars. Ha images of Hn 
regions ionized by OB associations in the Large 
Magellanic Cloud (LMC) have revealed large (50- 
150 pc) gas shells around the OB associations and 
numerous small (< 15 pc), ring- like nebulae out- 
side the large shell (see, e.g., Davies, Elliott, & 
Meaburn 1976). These small nebulae are candi- 
date wind-blown bubbles of single O or B stars 
(Weis & Duschl 1999), and offer perhaps the sim- 
plest test of the standard, adiabatic bubble model. 

We have investigated the kinematics and physi- 
cal nature of three ring-like nebulae in N44 (Hcnizc 
1956) in the LMC. N44 is a bright Hn complex 
with a 44 pcx67 pc superbubble at the center and 
several compact Hn regions along the shell rim. 
In the surrounding area, small dense H n regions 
as well as extended diffuse nebulae are also present 
at distances up to 160 pc. Figure 1 shows these 
features and their identifications given by Hcnizc 
(1956). Three OB associations exist in N44 (Lucke 
& Hodge 1970): LH 47 in the central shell, LH 48 
in the compact Hn region N44I, and LH 49 par- 
tially embedded in the bright Hn region N44D. 
The three nebulae that we have studied are N44F, 
N44J, and N44M. N44F is a bright, circular Hn 
region at the northwest rim of the superbubble and 
is adjacent to a bright filament connected to the 
main nebula. In contrast, N44J and N44M, ap- 
pearing respectively on the northern and eastern 
parts of N44, are fainter and relatively isolated. 
None of these three nebulae have been studied 



since their original identification by Henize (1956). 

We have obtained ground-based CCD images 
and high-dispersion cchellc spectra of N44F, 
N44J, and N44M to study the morphology and 
kinematics of these nebulae. Additionally, we have 
used Hubble Space Telescope (HST) WFPC2 im- 
ages of N44 F in the Ha and [S n] lines to analyze 
its ionization structure. In this paper, we first 
describe these observations in §2, then present the 
analysis and results for N44F in §3 and N44 J and 
N44M in §4. We present our conclusions in §5. 

2. Observations 

The datasets used in this study include: (1) 
CCD images taken with the Curtis Schmidt and 
0.9 m telescopes at Cerro Tololo Inter-American 
Observatory (CTIO), (2) CCD images taken with 
the HST WFPC2, and (3) high-dispersion cchclle 
spectra taken with the CTIO 4 m telescope. The 
high-dispersion spectra are useful to diagnose ex- 
panding shells that are not morphologically iden- 
tifiable. 

2.1. Ground-based CCD Images 

CCD images of N44 were taken with the Cur- 
tis Schmidt telescope at CTIO in 1993 Decem- 
ber. The detector was a front-illuminated Thom- 
son 1024x1028 CCD camera with 19 ^m pixels, 
or l'/835 pixel -1 . Images were obtained with 
Ha (A c =656lA, AA=26A), [Om] (A C =5010A, 
AA=50A), and [Sn] (A C =6720A, AA=50A) fil- 
ters for exposure times of 1800 s, 2700 s, and 
2700 s, respectively. The Ha filter includes not 
only the Ha line but also the neighboring [Nil] 
AA6548,6583 lines. At the LMC's radial velocity, 
-300 km s"\ the [Nil] A6583 line is red-shifted 
further away from the red edge of filter, but the 
[N n] A6548 line is red-shifted further into the Ha 
bandpass. Using the [Nil] /Ha intensity ratio mea- 
sured from our echellc spectra (see Section 2.3), 
we estimate that the [Nil] A6548 line emission 
contributes 2.4%, 4.3% and 3.3% of the total flux 
of N44F, N44J and N44M, respectively. This 
[N n] contamination has been corrected in our Ha 
flux measurements. 

Higher-resolution Ha images were taken with 
the 0.9 m telescope at CTIO in 2000 Decem- 
ber. The detector was a thinned, AR coated 
SITe 2048x2048 CCD with 24 ^im pixels, or 
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0"4 pixel- 1 . The Ha filter (A C =6559A, AA=64A) 
used transmits Ha and both [Nil] lines, therefore 
the "Ha" images are really Ha+[Nli] images, and 
thus are only used for morphology analysis. The 
angular resolution, determined by the FWHM of 
the stellar images in the field of view, was ap- 
proximately l'/3. To map the entire N44 complex, 
five fields were observed with exposure times of 
2x600 s for each field. 

All CCD images were bias-subtracted and flat- 
fielded in a standard way using IRAF software, 
and multiple frames were combined to remove 
cosmic rays and improve the S/N ratio. Fig- 
ure 1 shows the Ha mosaic of the N44 complex, 
and Fig. 2 presents close-ups of N44F, N44 J and 
N44M in the Ha, [Om], and [Sn] lines. 

2.2. HST WFPC2 Images 

HST WFPC2 images of N44F were acquired 
via the Hubble Observation Problem Report 7 
(HOPR) observations of the Cycle 6 program 
6698. In these observations, made with the 
F656iV filter for 2x500 s and the F673A filter for 
2x600 s, N44F was included in the Wide-Field 
Camera 3 (WF3). The F656N filter, centered at 
6563.7 A with a FWHM of 21.4 A, includes the 
Ha line and the neighboring [Nil] lines. Using 
the [Nil] /Ha intensity ratio measured from our 
echelle spectra (see Section 2.3), we estimate that 
the [Nil] A6548 line emission contributes 2.4% of 
the flux in N44F. This [Nil] contamination has 
been corrected in our Ha flux measurements. The 
F673iV filter, centered at 6732.2 A with a FWHM 
of 47.2 A includes the nebular lines of [Sn] AA 
6716,6731. 

The calibrated WFPC2 images were produced 
by the standard HST pipeline. We processed them 
further with IRAF and STSDAS routines. The 
images taken with the same filter were combined 
to remove cosmic rays and to produce a total- 
exposure map. The combined [S n] and Ha images 
were then each corrected for the intensity- and 
position-dependent charge transfer efficiency as 

7 The original observations of N44 were made with an incor- 
rect orientation, resulting in a loss of 1/3 of the field of 
view. To compensate the loss, HOPR observations were 
granted. To facilitate the scheduling of these HOPR ob- 
servations and to maximize the scientific yield, a different 
pointing center in N44 suitable for arbitrary orientations 
was chosen. 



explained by Holtzman et al. (1995). Following the 
recommended procedures for narrowband WFPC2 
photometry 8 , we find the rectangular widths of the 
,F656iV and F673A filters to be 28.3 A and 63.3 
A, respectively, and use these widths to convert 
the flux densities into fluxes. 

Figure 3 shows the WFPC2 images of N44F. To 
show the excitation variations, the [Sn]/Ha ratio 
map of N44 F is also presented, in addition to the 
Ha and [Sn] images. 

2.3. CTIO 4 m Echelle Spectra 

Echelle spectra were obtained with the CTIO 4 
m telescope in 2000 December. The spectrograph 
was used with a 79 line mm -1 echelle grating and 
the long-focus red camera. We observed a sin- 
gle order centered on the Ha line by inserting a 
post-slit Ha interference filter and replacing the 
cross-dispersing grating with a flat mirror. A Tek 
2048x2048 CCD with 24 fim pixel" 1 was used to 
record the images. This provided a spectral sam- 
pling of ^0.08 A pixel" 1 and a spatial sampling 
of 0'.'26 pixcr 1 . The wavelength coverage, limited 
by the Ha filter and the echelle order, was 125 
A; the spatial coverage, limited by the optics, was 
^200". The angular resolution, determined by the 
FWHM of the seeing, was approximately 1". A 
slitwidth of 1'.'65 was used. The resultant spectral 
resolution, measured from the Th-Ar lamp lines, 
was about 13 km s" 1 FWHM. 

Each of the three nebulae, N44F, N44J, and 
N44M, was observed with an east-west oriented 
slit. The exposure times were 600 s for N44F, but 
only 300 s for N44 J and N44M. The spectral lines 
detected include nebular Ha and [Nil] lines, and 
telluric Ha and OH lines (Osterbrock et al. 1996). 
We have adopted rest wavelengths of 6562.7885 
A, and 6548.080 A and 6583.454 A for the Ha 
and [Nil] AA6548,6583 nebular lines, respectively 
(Spyromilio 1995). 

3. Physical Structure of N44 F 

We will first describe the overall morphology of 
N44F, examine its excitation and ionization struc- 
ture, and derive its rms electron density. Special 
emphasis is given to the prominent dust "pillars" 
where star formation might be taking place. Next, 

available at http://www.stsci.edu/instruments/wfpc2/Wfpc2_faq/wfpc2_nrw_phot. 
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we describe the kinematic features detected in the 
echelle spectra, fit Gaussian components to the 
line profiles, and use the resultant velocities to 
identify expanding structures and to determine ex- 
pansion velocities. Finally, we combine the mor- 
phological and kinematic information and present 
an integrated view of this H n region. Throughout 
this paper, we will adopt a distance of 50 kpc to 
the LMC (Feast 1999). 

3.1. Morphology and Density 

N44F is a bright circular nebula nestled on a 
bright filament at the northwest rim of the su- 
pcrbubblc in N44. In low-resolution Ha images, 
N44F appears to have a double- ring structure 
with a very bright region in the northwest. The 
details of these features are resolved in the HST 
WFPC2 images. The outer ring, of radius ~ 22", 
is surrounded by dusty features. Finger- or pillar- 
like dust columns, similar to those in the Eagle 
Nebula (=M16; Hester et al. 1996), are also seen 
at the periphery, with two especially prominent 
ones at the west, close to the brightest region of 
N44F. These dust columns are revealed in detail 
by the HST images, and will be discussed more 
extensively in Section 3.2. The inner ring, of ra- 
dius ~ 6", is resolved into diffuse emission with- 
out sharp features. The western part of this ring 
merges into the brightest region of N44 F. A group 
of bright stars are located within the central ring, 
but they are substantially offset towards the bright 
emission region at the northwest. Unfortunately, 
none of the stars in the vicinity of N44 F have been 
studied photometrically or spectroscopically. 

The [S n] images of N44 F show a very differ- 
ent morphology: only the outer ring is visible and 
it is mostly composed of sharp filaments. The 
pillar- like features are also present and stars at the 
tips of two of them are clearly seen (see Fig. 4). 
The [Sn]/Ha ratio in the central region is very 
low, about 0.1. The outer ring possesses a higher 
[Sn]/Ha ratio, with values ranging from 0.2 to 
0.4. Enhanced [Sn]/Ha ratios are often associ- 
ated with shell structures or high velocity shocks. 
As we do not find the latter in the nebula (see Sec- 
tion 3.3), the somewhat elevated [Sn]/Ha ratio is 
likely to result from the low ionization parameter 9 

'The ionizing flux derived further can be used to estimate 
the ionization parameter U oc (Q(H°) n Z 2 ) 1 / 3 , where 



caused by the shell morphology. The values of 
[Sii]/Ha are in fact similar to those seen in other 
shell structures (e.g., Lasker 1977; Hunter 1994). 

Using the Curtis Schmidt images, we also exam- 
ine the [O in] morphology of N44 F. Because of the 
poorer resolution, only the overall morphology can 
be described. The [O in] image of N44F looks very 
similar to that in Ha, except for the bright Ha 
outer ring, which disappears almost completely in 
[Oiii]. We have used the HST WFPC2 observa- 
tions of N44C (Garnett, Galarza, & Chu 2000) 
to calibrate the observed [O in] /Ha ratio derived 
with the Curtis Schmidt images. The [Oni]/Ha 
ratio of N44F is relatively high, about 0.5-1.2 
inside the ring, with the highest value near the 
star(s). The outer ring has the lowest [O in] /Ha 
ratio, about 0.25. The integrated [Om]/Ha ratio 
for the whole nebula is about 0.55. The high ratio 
near the central region indicates that the ionizing 
source must be an early-type massive star. 

The WFPC2 Ha image can be further ana- 
lyzed to estimate the rms density of N44F. A 
color excess of E(B— V)=0.11 mag for N44 (Oey 
& Massey 1995) is adopted in the extinction cor- 
rection. The emission measure (EM) derived from 
the extinction-corrected surface brightness (Pcim- 
bert, Rayo, & Torres-Peimbert 1975) ranges from 
7xl0 3 cm~ 6 pc in the faintest interior regions of 
N44F to 4xl0 4 cm" 6 pc at the tips of the "pil- 
lars" . To derive the density from the EM, we can 
assume two simple geometries. First, the density 
can be estimated by assuming the depth to be 
equal to the diameter of the nebula: for the cen- 
tral region, we find a value ranging from 25 cm~ 3 
to 35 cm~ 3 . Second, N44F is limb-brightened 
and we can therefore estimate the density by as- 
suming a constant density spherical shell with ra- 
dius R and thickness AR. This geometry predicts 
a unique radial surface brightness profile for any 
AR/R based on the varying path length through 
the shell material. Comparing these simple mod- 
eled profiles to the observed Ha surface brightness 
profiles near the shell rim, we find a shell thickness 
AR/R of ^0.25. Using the peak surface bright- 
ness and the length of the longest line of sight, 

Q(H°) is the ionizing flux in photons s -1 , n the density 
in cm" 3 , and / the filling factor. Using Q(H {, )= l.OxlO 49 

photons s -1 , n=40 cm -3 and /=1 — ^1 — ^-*J ~0.6, we 

find a value of ~ 10~ 3 for the ionization parameter. 
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2y/2RAR — (Ai?) 2 , we then derived a density in 
the shell of 40 cm~ 3 . On average, the diffuse gas 
of N44 F probably possesses a density of a few tens 
cm~ 3 . 

Using the Curtis Schmidt Ha image and the 
flux calibration from Kennicutt & Hodge (1986), 
we measured over the entire N44 F nebula an Ha 
flux of 3.6xl0~ n erg cm~ 2 s _1 . For a distance 
of 50 kpc, the Ha luminosity of N44F is then 
1.4xl0 37 erg s -1 , after correcting for an extinc- 
tion of E(B— V)=0.11 mag. This Ha luminos- 
ity requires an ionizing flux of Q(H°) — 1.0x 10 49 
photons s _1 . A second evaluation of the extinc- 
tion correction has been made by comparing radio 
data and the Curtis Schmidt Ha image: a mean 
value E(B-V) of 0.28 mag for N44F, N44 J, and 
N44M is found (J. R. Dickel, private communica- 
tion), and this results in a Q(H°) of 1.4x 10 49 pho- 
tons s _1 for N44F. This value suggests that the 
total ionizing power of the star(s) within N44F is 
equivalent to a single 06.5 V-07.5 V star (Pana- 
gia 1973; Vacca, Garmany, & Shull 1996; Schaerer 
& de Koter 1997). This is compatible with the 
rather high [O in] /Ha ratio: the photoionization 
models of Stasihska & Schaerer (1997) for CoStar 
stellar atmosphere models suggest that the ob- 
served integrated [O in] /Ha ratio of N44F could 
be produced by the equivalent of an 08 V-09 
V star 10 . The comparison to the photoionization 
models should however be interpreted with cau- 
tion because the models do not assume a shell 
structure, which affects the excitation and ioniza- 
tion structure. The dominant spectral type could 
therefore be somewhat earlier than inferred from 
the modeled excitation. 

3.2. Pillar-like Dust Columns 

The dark clouds of N44 F exhibit new examples 
of the intriguing finger- or pillar-like structures, 
similar to those reported in the Eagle Nebula (Hes- 
ter ct al. 1996). The most prominent two are em- 
bedded in the bright northwest region (see Figs. 3 
and 4). Several additional features of this type 
are present in the south and possibly in the south- 
west, but they are fainter. In Fig. 4, we present 

'We have compared N44F, N44J, & N44M to the closest 
model available in Stasiriska & Schaerer (1997), i.e. a 
spherical nebula of density 10 cm -3 and metallicity 0.25 
Z©, ionized by one single star. 



the most striking of these structures in the Ha and 
[S n] images as well as in the ratio map. The dust 
columns in N44F arc 5" to 10" long and 1" to 2 
" wide, corresponding to 1 to 2 pc long and 0.25 
to 0.5 pc wide at the LMC's distance. These sizes 
and morphologies are similar to those seen in 30 
Doradus (Scowen et al. 1998) and M16 (Hester et 
al. 1996). Two of these columns harbor stars at 
their tips, as indicated by arrows in Fig. 4. Hes- 
ter et al. (1996) suggested that stars at the tips of 
pillars could be young stellar objects which were 
formed in the molecular cloud and finally emerge 
after the photoevaporation of most of their parent 
cloud. However, the Hester ct al. (1996) model of 
pre-existing condensations in the molecular cloud 
does not provide the unique conditions for pro- 
ducing photoionized pillars. Hydrodynamical sim- 
ulations by Williams, Ward-Thompson, & Whit- 
worth (2001) illustrate that these structures can 
be formed by a variety of initial conditions in the 
structure with inhomogeneities in density and/or 
radiation field, but current observations do not yet 
permit discrimination among these different mod- 
els. It is therefore of great interest to identify ad- 
ditional examples of such objects. 

Other important characteristics of these dust 
columns are their ionization structure and their 
density. All of the dust columns point towards 
the central ionizing source, i.e., a bright star or 
a small cluster between the center and the north- 
west edge of N44F. As can be seen in Fig. 5, the 
[Sn] peak is separated from the Ha peak by ap- 
proximately 071, i.e., 0.024 pc, with Ha peaking 
slightly closer to the source of UV radiation which 
is photoevaporating the cloud. The shape of these 
profiles is very similar to those observed in the 
pillars of 30 Doradus, or in ground-based images 
of M16 (Scowen et al. 1998). Unfortunately, no 
WFPC2 [O in] images of N44 F are available, thus 
we can not undertake a complete analysis using 
photoionization models or constrain further the 
properties of the ionizing source. Following the 
method used in section 3.1, wc estimate an rms 
density larger than 300 cm -3 at the tip of one of 
the pillars. This density is compatible with that 
determined by Scowen et al. (1998) in 30 Doradus, 
700 cm~ 3 . These pillars are thus photoevaporat- 
ing clouds, whose evaporated matter enhances the 
density and brightness of the diffuse material seen 
in the northwest part of the nebula. 
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3.3. Kinematic Properties 

The kinematic properties of N44 F are derived 
from the high-dispersion cchclle observations of 
the Ha and [N n] A6583 lines. The [N n] line has 
a smaller thermal width than the Ha line, thus it 
shows the line splitting better although it is nois- 
ier. We have extracted Ha and [Nil] line profiles 
and fitted them with Gaussian components. The 
resultant velocity-position plot (Fig. 6) shows that 
N44F is an expanding shell with a radius of 22", 
coincident with the outer ring described in §3.1. 
Both lines indicate a shell expansion velocity of 
^12 km s _1 . An example of the [Nil] and Ha 
profiles are presented in Fig. 7. The approach- 
ing (blueshifted) component is fainter, suggesting 
a lower density. Further east, outside N44F, an 
additionnal blue component is present, but it is 
probably associated with the complex filamentary 
structure at the outskirts of N44 and is not directly 
related to N44F. 

3.4. Test of the Standard Bubble Model 

N44F has a classical ring-nebula morphology, 
and its kinematics confirms an expanding shell, 
most likely created by the wind of the early-type 
massive star(s) within this bubble. We can now 
test whether the observed expansion dynamics of 
this bubble are consistent with those expected 
from the standard, adiabatic bubble model (e.g., 
Weaver et al. 1977). With the observed radius 
R in units of pc and expansion velocity V in 
units of km s _1 , we can compute the dynamical 
timcscale t$ in units of 10 6 yr using t§ = Q.6R/V, 
and further determine the ratio of stellar wind 
luminosity L 36 in units of 10 36 ergs s _1 to am- 
bient density no in units of H-atom cm~ 3 , using 

For a radius of 5.3 pc and an expansion ve- 
locity of ~12 km s _1 , N44F has a dynamical 
timescale of 2.7xl0 5 yr and a L 36 /n ratio of 0.016 
ergs s _1 cm 3 . If we assume that the shell consists 
of mass swept up from the volume, the ambient 
density, no, can be estimated by n = 3n-^, 
where n is the rms density of the shell, AR the 
shell thickness, and R the shell radius. For n ~40 
cm~ 3 and ^ = 0.25 (see Section 3.1), we derive 
an ambient density of ~30 cm~ 3 and a wind lumi- 
nosity of ^5xl0 35 ergs s _1 . This wind luminosity 
is compatible with that expected from the spectral 



type determined in Section 3.1, i.e., an 07V star 
(de Jager et al. 1988; Prinja, Barlow, & Howarth 
1990). One uncertainty in these estimates is the 
shell thickness; however, even if the shell thickness 
is off by a factor of two, the predicted wind lumi- 
nosity will still be compatible, within the errors, 
with an 07V star. 

Thus, the dynamics of N44F appear to be con- 
sistent with the standard adiabatic model for bub- 
ble evolution, within uncertainties generated by 
the complex interstellar environment. This con- 
trasts with the results from superbubbles, which 
show evidence of retarded growth (e.g., Oey 1996). 
To investigate further this discrepancy, it is worth 
considering an alternative bubble model. 

The presence of dust pillars around N44F 
clearly indicates a highly inhomogeneous ISM. 
The photoevaporation of the dense clumps may 
serve as an energy sink and govern the dynamic 
evolution of this H n region. Since photoevapora- 
tion processes take place in the early evolutionary 
stages of a bubble/superbubble, the energy loss by 
photoevaporation may explain the smaller-than- 
expected sizes of previously studied superbubbles 
and Wolf-Rayet ring nebulae (see §1). Therefore, 
we will now consider models of bubbles formed in a 
cloudy medium by McKee, Van Buren, & Lazareff 
(1984). In their models, only bubbles blown by 
stars with weak winds will follow the Weaver et 
al. model, while stars with stronger winds will 
produce bubbles whose expansion is governed by 
the dynamics of photoevaporating clouds. They 
define a dimcnsionless wind luminosity L* w = 
L w /(1.26 x 10 36 (5| 9 /n m ) 1 / 3 ) and a homogeniza- 

1/7 

tion radius R h (pc) = 0.009 i 4/7 (ff^ , where 

L\y is the wind luminosity in erg s , S49 the 
ionizing flux in photons s _1 , n m the mean cloud 
density in cm~ 3 and t the dynamical timescale in 
years. Strong winds possess > 1 and blow 
bubbles of radius ~ Rh- 

Using the observed parameters of N44 F (with 
n m =ri{)) and assuming a stellar wind luminos- 
ity for an 07V star, we estimate ~1 and 
Rh ^4.5 pc, close to the observed bubble's radius. 
This seems to suggest that the dynamics of photo- 
evaporation actually governs the bubble's expan- 
sion. The agreement of the data with both models 
could be a simple coincidence, or maybe the bub- 
ble has just encountered the cloudy medium and 
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begun the photoevaporation processes, in which 
case there may not be a large deviation from the 
Weaver et al. model. Note, however, that all these 
conclusions hinge upon the validity of the stellar 
wind luminosity, which is based on an assumed 
spectral type of the central star. Future study of 
the actual stellar content of N44F is needed to 
place stronger constraints on the applicable bub- 
ble models. 

4. Physical Structure of N44 J and N44 M 

4.1. N44J 

Situated in the northernmost part of the N44 
complex, N44J is a slightly elongated, ring-like 
nebula with a size of 37"x32"(or 9.0 pcx7.8 pc). 
The 0.9 m Ha image shows a bright rim around 
dusty features at the northwest boundary. A 
bright central star is visible in the Ha and [S n] 
images (see Fig. 2). In Ha and [Sn], the nebula 
exhibits modest limb-brightening, with the edge- 
to-center brightness ratio <2, while its central part 
displays a rather uniform surface brightness. The 
[Sn]/Ha ratio is very similar to that of N44F, 
ranging from 0.25 to 0.4. In contrast, the [Om] 
image of N44 J shows emission only near the cen- 
tral star. The [Om]/Ha ratio is 0.4 near the star, 
and only 0.05-0.08 elsewhere, significantly lower 
than the corresponding value in N44F. The inte- 
grated [O III] /Ha ratio of N44 J is about 0.12. 

The observed Ha flux of N44J is 5.7xl0~ 12 
erg cm~ 2 s _1 . This requires an ionizing flux 
Q(H°) of 1.6xl0 48 photons s -1 when adopting 
E(B-V)=0.11 mag, or 2.3xl0 48 photons s" 1 for 
E(B— V)=0.28 mag. These values are expected for 
an O9V-B0V type ionizing star (Panagia 1973; 
Vacca, Garmany, & Shull 1996; Schaerer & de 
Koter 1997), which is fully consistent with that 
suggested by the [O in] /Ha ratio of the nebula 
(Stasihska & Schaerer 1997). Assuming a spher- 
ical geometry, we derive an rms electron density 
of ~25 cm~ 3 , slightly less dense than N44F. The 
ionized gas mass of N44 J is then ^200 M Q . 

The echelle spectrum of N44J shows a kine- 
matically quiescent Hn region. The centroids 
of the Ha and [Nil] lines are centered at 292 
km s _1 , with an rms variation of 3 km s _1 (see 
Figs. 6 and 7). No line splitting is seen. The ob- 
served FWHMs are 27±1 km s" 1 for Ha and 16±2 
km s _1 for [Nil]. These observed line widths im- 



ply a turbulent FWHM of about 10 km s _1 for 
Ha and 8 km s _1 for [Nil] if we take into account 
the instrumental profile and the thermal broad- 
ening at 10 4 K. The turbulence can thus explain 
the small variations of the centroid velocity. The 
difference between the turbulence widths is prob- 
ably caused by a difference in the emission regions 
of the two lines. If N44J is actually expanding, 
the expansion velocity should be very small, <8 
km s _1 . For such small expansion velocity, no 
strong compression is expected, hence no sharp 
filaments should be seen. 

Finally, we may compare the observations with 
the standard, adiabatic bubble model. As the den- 
sity in the ring-like nebula is a few tens cm~ 3 , 
we will assume in our calculations a density of 10 
cm~ 3 in the ISM. For an expected wind luminos- 
ity of ~3xl0 34 ergs -1 for an 09.5V star (de Jager 
et al. 1988; Prinja, Barlow, & Howarth 1990) and 
with the observed 17'.' 5 radius of N44 J, the Weaver 
et al. (1977) model predicts an expansion velocity 
of about 8 km s -1 for N44J. This expected ex- 
pansion velocity is just at the limit implied by the 
velocity widths of the Ha and [Nil] lines. Rais- 
ing the ambient density by a factor of 2 can lower 
the expected expansion velocity to 6.5 km s _1 , 
which is more consistent with the observed limit. 
We conclude that the observed bubble dynamics of 
N44 J does not greatly contradict the model pre- 
dictions. 

4.2. N44M 

On the southeast part of N44, another small 
ring-like nebula is present: N44 M. It is the faintest 
nebula of our sample. It shows a ring with a ra- 
dius of ~22" (or 5.3 pc), but its north bound- 
ary is not well defined (see Fig. 2). The sharp 
southwest boundary suggests the presence of dust 
clouds. Several bright stars appear in its interior. 
N44 M possesses an excitation structure very sim- 
ilar to N44 J. The [Sn]/Ha ratio of N44M is com- 
parable to that in N44F, but the [O III] /Ha ratio 
of N44M is much lower, ranging from 0.07 to 0.16 
in the nebula and rising to 0.4 near the stars. The 
integrated [O III] /Ha ratio of N44 M is about 0.18. 

With an observed Ha flux of 5.3xl0~ 12 erg cm~ 2 s _1 
and the two evaluations of the extinction, the re- 
quired ionizing flux Q(H°) is 1.5-2. lxlO 48 pho- 
tons s _1 . This suggests an 09-BOV type ionizing 
star (Panagia 1973; Vacca, Garmany, & Shull 
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1996; Schaerer & de Koter 1997), which is again 
compatible with the [O in] /Ha ratio (Stasihska & 
Schaerer 1997). If we assume a spherical geom- 
etry, we hnd that the rms density is about ~17 
cm~ 3 and the estimated mass of the nebula ~275 
M . 

The echelle spectrum of N44 M shows only one 
component, centered at 298 km s _1 with an rms 
variation of 1 km s _1 , and a FWHM of 28±1 
km s _1 in the Ha line and 20±4 km s _1 in the 
[Nil] line (see Figs. 6 and 7). The turbulent 
FWHM arc then 12 km s" 1 for Ha and 14 km s" 1 
for [N n] . With an observed radius of ~22", a den- 
sity of 10 cm~ 3 for the ISM and a 09.5V ionizing 
star (see Section 4.1.), the standard bubble model 
predicts an expansion velocity of ~7 km s _1 for 
N44M, which appears consistent with the upper 
limit of the echelle observations. 

The kinematic characteristics of N44J and 
N44M are very similar, despite the somewhat 
different morphology. Neither Hn region con- 
tains shells expanding at velocities greater than 
~8 km s _1 . In the future, additional observations 
may enable a better comparison of these nebu- 
lae with the standard model. Precise photometry 
and/or stellar spectroscopy would give accurate 
informations about the stellar content of these 
nebulae, while nebular spectra with higher reso- 
lution should more strongly constrain the kine- 
matics of N44 J and N44M. Thus, improved data 
is needed to more stringently test the standard 
model for N44 J and N44M. 

5. Conclusion 

Massive stars inject a large amount of energy 
into the ISM through their ionizing radiation, fast 
stellar winds, and SNe. The stellar winds, which 
dominate the early mechanical feedback, are able 
to sweep the ambient medium into an expanding 
shell and create a cavity around the parent stars. 
H II regions ionized by massive stars are therefore 
thought to evolve from compact Hn regions to 
ring-like nebulae. 

We have studied three candidate wind-blown 
bubbles in the N44 complex: N44F, N44J, and 
N44M. They all have diameters of ~10 pc. Their 
ionizing fluxes and [O in] /Ha ratios suggest that 
their ionizing stars have spectral types of roughly 
07V, 09.5V and 09.5V, respectively. N44F has 



an expansion velocity of 12 km s _1 , while the 
other two have expansion velocities < 8 km s _1 . 
The stellar wind luminosity implied by the ionizing 
stars and the observed bubble dynamics appear 
consistent with that expected in the standard, adi- 
abatic bubble model within the limits of observa- 
tional uncertainties. Similar results were found by 
Oey & Massey (1994) for two bubbles around sin- 
gle star in M33. This is in sharp contrast to the 
anomalous kinematics seen in superbubbles and 
Wolf-Rayet bubbles. We caution, however, that 
our observed expansion velocities may not pro- 
vide the most stringent constraints when compar- 
ing with models, because the expected expansion 
velocities of N44 J and N44 M are comparable to 
or smaller than the convolution of the isothermal 
sound velocity (~10 km s -1 ) and the instrumental 
FWHM (13 km s _1 ). Deeper, higher-resolution 
spectroscopic observations of the [N n] or [O in] 
lines are needed to pin down the expansion veloc- 
ity for a more definitive test of models. 

We also report the discovery of several dust pil- 
lars in N44F, which are similar to those found in 
the Eagle Nebula (Hester et al. 1996) and 30 Do- 
radus (Scowen ct al. 1998). These structures are 
indicative of the evaporation of molecular clouds 
by the ionizing source. The examples in N44F are 
especially well-suited to investigating the process 
and role of photoevaporation in star formation and 
shell evolution, because of their location within 
such a well-defined, wind-blown bubble. We found 
that photoevaporation, as evidenced by the pillar 
structure, may be dynamically important for the 
shell evolution in N44 F, slowing down the expand- 
ing shell and inhibiting its growth. These pho- 
toevaporation processes may thus undermine the 
seeming agreement between the observed bubble 
dynamics and the simple adiabatic bubble model. 
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Fig. 1.— CTIO 0.9 m Ha (+[Nn]) CCD im- 
age of N44. The nebular components defined by 
Henize (1956) and the OB associations cataloged 
by Lucke & Hodge (1970) are marked. 



Fig. 2. — Ground-based images of N44F (left), 
N44J (middle), and N44M (right). The first 
row contains Ha (+[Nn]) images taken with the 
CTIO 0.9 m telescope. The second, third, and 
fourth rows are Curtis Schmidt images taken in 
Ha (+[Nn]), [Oin], and [Sn] filters, respectively. 
The field-of-view (FOV) of each image is 92"x92". 
North is up and east to the left. 



Fig. 3.— HST WFPC2 Ha (+[Nn]) image (left), 
[Sn] image (middle) and [Sn]/Ha ratio map 
(right) of N44F. The rectangular boxes show the 
position of the close-up images presented in Fig. 
4. The bars below the images show the greyscales 
for surface brightness (in units of erg cm~ 2 s _1 
arcsec~ 2 ) or line ratio. North is up and east to 
the left. 



Fig. 6. — Top panel: Echellogram of the [Nil] 
A6583 line for N44F. Bottom panels: Radial 
velocity-position plots of the echelle slit positions 
in N44F, N44 J, and N44M for both Ha and [Nil] 
lines. The radial velocities are heliocentric. The 
Ha velocity components are plotted in filled sym- 
bols, and the [N II] in open ones. The position 
axis is defined with the center of the nebula at 
the origin (0") and increasing towards west. The 
dot-dashed line in each plot indicates the range 
covered by the nebula. Arrows show the positions 
of the profiles presented in Fig. 7. 



Fig. 4. — Close-up images of dust columns seen in 
N44F through Ha, [Sn], and [Sn]/Ha (for exact 
position, see Fig. 3). Left: dust columns situated 
at the west rim of N44F (FOV=18"xl8"). Right: 
dust columns situated at the south rim of N44F 
(FOV=14"xl4"). The arrows point to the stars 
at the tips of the pillar-like dust columns. 



Fig. 5. — Left: Ha image of the "pillars" situated 
at the west part of N44F. The rectangular box 
and the arrow show the position and direction of 
the cut presented at the right. Right: Ionization 
structure along the rectangular aperture. [S n] and 
Ha intensities arc normalized to the peak, while 
the [Sn]/Ha ratio is the actual value. Note that 
the Ha and [S n] peaks are separated by one pixel, 
i.e. O'.T. 



Fig. 7. — Ha and [Nil] profiles at the positions 
indicated in Fig. 6. 
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